High-field vortices in Josephson junctions with alternating critical current density 
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We study long Josephson junctions with the critical current density alternating along the junction. 
New equilibrium states, which we call the field synchronized or FS states, are shown to exist if the 
applied field is from narrow intervals centered around equidistant series of resonant fields, Hm- 
The values of Hm are much higher than the flux penetration field, Hs- The flux per period of the 
alternating critical current density, ^i, is fixed for each of the FS states. In the m-th FS state the 
value of is equal to an integer amount of flux quanta, 0^ = m(f)Q. Two types of single Josephson 
vortices carrying fluxes 00 or/and (j}o/2 can exist in the FS states. Specific stepwise resonances in 
the current- voltage characteristics are caused by periodic motion of these vortices between the edges 
of the junction. 

PACS numbers: 74.50.+r, 74.78.Bz, 74.81.Fa 

Keywords: vr -shifted Josephson junction, Josephson vortex 
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Josephson tunnel structures arranged in a sequence of 
interchanging - and tt - shifted fragments (see Fig. TO) 
are a subject of growing interest [J [3, H, Q, d, 0, H, 
i, 1 113, 111 E, 111, [3 111- Currently, the proper- 
ties of these complex systems are thoroughly treated for 
superconductor-ferromagnet-superconductor (SFS) het- 
erostructures [ij, |2|, (sj, |3, [1| and asymmetric grain bound- 
aries in thin films of high-T^ superconductor YBCO 
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It has been predicted that in the SFS Josephson junc- 
tions the TT-shift in the phase difference ip between the 
superconducting banks is caused by the ferromagnet in- 
terlayer [J, [|. This prediction was confirmed in re- 
cent experiments with SFS tt- shifted junctions and SFS 
heterostructures of interchanging - and tt - shifted frag- 
ments [1, 13, @- The asymmetric grain boundaries in 
YBCO thin films are arranged in series of facets with vari- 
ety of orientations and lengths Z ~ 10 — 100 nm [l3| . This 
spatial structure in conjunction with the d-wave sym- 
metry of the order parameter results in grain boundary 
junctions with interchanging - and tt - shifted fragments 
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The critical current density, jc{x), changes sign at each 
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contact between 0- and tt- shifted fragments {x axis is 
along the junction). This results in a variety of dra- 
matic changes in Josephson properties, in particular, the 
field dependence of the maximum supercurrent across the 
junction, 7^, is strongly affected by the alternating jc{x) 
@, 0, [l3, EE]- First, the value Im(,Ha) is significantly 
suppressed at low fields, \Ha\ Hi = 0o/2AZ, where A 
is the London penetration depth and I is the period of 
jc{x). Second, unlike the standard Fraunhofer pattern 
with a major peak at Ha = 0, two major side-peaks are 
observed at high fields. Ha = ±7Ji, where Hi ^ Hg and 
Hs is the field of first flux penetration [7|. 

In this letter we find a series of new equilibrium field 
synchronized (FS) states existing if the applied field is 
from narrow intervals AHa ^ Hs centered at the reso- 
nant fields H„i = ±miJi, where m ^ is an integer. It 
is shown that in the m-th FS state the flux per period 
of the alternating critical current density, 4>i, is fixed and 
is equal to an integer number of flux quanta, (pi — m(j)o . 
We demonstrated that two types of high-fleld {Ha ^ Hg) 
Josephson vortices with fluxes 0o or/and (f)o/2 can exist 
in the FS states. We find that periodic motion of these 
single high-field vortices cause stepwise resonances in the 
IV curves similar to the zero-field resonances in Joseph- 
son junctions of conventional superconductors (l6l |. 

We begin with a qualitative treatment of the FS states 
using one harmonics model for the tunneling current den- 
sity j = jc{x) svaip, where jc{x) = ji sin(27ra;/Z), L = Nl, 
L is the length of the junction, and ^ 1 is an inte- 
ger. Assume that the junction is in one of the FS states 
and flux (pi is fixed. Since we have many vortices in the 
junction {N • m 3> 1), the field is almost uniform and the 
phase (p{x) takes the form 



FIG. 1: Schematic diagrams of: (a) tunnel junction arranged 
in series of interchanging - and tt -shifted fragments; (b) spa- 
tial distribution of alternating critical current density. 



(/?(x) = 27r^y +i'{x).. 
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where the phase ^{x) is a smooth function with the typ- 
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ical length scale ^ I and ^ 1- Then we have: 



j{x) ^ ji sin 
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In general, this current density alternates rapidly with a 
typical length scale < In this case the coarse-grained 
approach is the right tool to describe the smooth phase 
ip^x) [l3|- If we average j(x) over a distance I, then 
the coarse-grained tunneling current density, j^, is zero. 
This is indeed true in all cases but one. If the junction 
is in the FS state with 1/)^ = ±0Oj then coarse-graining of 
Eq. ([2]) leads to a nonzero result 



= 0.5 ji ain{ip ± 7r/2). 



(3) 



It is worth mentioning that ~ ji and the final form 
of the dependence of on the smooth phase ip coincides 
with a 7r/2 -shifted current-phase relation of a Josephson 
junction of conventional superconductors. 

Similar calculation of the coarse-grained current den- 
sity in the low-field region (Ha <€. Hg) leads to 
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sin 2ip < ji , 



,47rAi 

where the Josephson length is given by 



Ai = v/c0o/167r2Aji > I . 



(4) 



(5) 



Comparing Eqs. ([2]), ([3]), and (g]) we find that if the 
phase factor, sin ip, modulations caused by the field and 
the critical current density modulations caused by the 
structure of the junction are synchronized — m(/)o), 
then the critical current density jc{x) is significantly en- 
hanced (see Fig. [2]). The widths of the intervals AHa ~ 
Hs ■€i Hi of existence of the FS states are defined by 
minimization of the free energy of the junction. 

Consider now a Josephson junction with A <C ^ ^ Aj, 
where Aj is the local Josephson penetration depth 



Aj = y/c^o/l67:^X{\jc\), (6) 
and the averaging over the junction length is defined as 



(/)- / dxf{x)/L. 



(7) 



In this case the static spatial distribution of the phase 
ip{x) is given by 



Ajiyj" — ic{x) sin ip = , 



(8) 



where ic{x) — jc{x)/ {\jc\) is the dimensionless critical 
current density. Next, we expend ic{x) in Fourier series 
and obtain instead of Eq. ([S]): 



sin = , 



(9) 
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FIG. 2: (a) Coarse-grained tunneling current density depen- 
dence on the applied field for the case of a stepwise critical 
current density, A*' — 40, and (p on x. Schematic diagrams of 
spatial distributions of tunneling current density for junctions 
that are: (bl) not in the FS states (small current loops tend 
to cancel each other); (b2) in the FS states (field synchronized 
local currents flow in the same direction). 



where i„ are the Fourier transforms of the function iJx) . 
Since I Xj, the phase (p can be written as 13l. Il7| 



If — 2'K(f)ix/ (j)ol -I- ^"(2;) + ^{x) 



(10) 



where ^j{x) is a smooth function with the typical length 
scale Aj ^ Z and |V'(2;)| ^ 1, £,{x) is a rapidly alter- 
nating function with the length scale ~ L In addition, 
we assume that (f(x)) = 0, and (|f(a;)|) ^ 1. Follow- 
ing 1^, ill we average Eq. © over the junction length 
and obtain the equation describing the smooth (coarse- 
grained) phase ip{x) in the m-th FS state 



Aj?A" = |i„i| sin(V' 
where Om = arg(i 
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Im) ~ Im sin 2(11} - ajn), (11) 
and 7,„ are defined by 
I 
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In general, complexity of Eq. pT|) can be significantly 
reduced. Indeed, one can estimate 7^ ~ {1/2ttXj)^ ^ 1, 
which means that "fm <C As a result the second term 
in the RHS of Eq. HT]) can be neglected and Eq. (fTTI) 
yields the sine-Gordon equation 



Ali," = sin(V' - 0„O . 



(13) 



where A^ = A,//-\/|i.,„|. The boundary conditions to Eq. 
([T3|) are given by the field at .t = and x ^ L: 
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It is worth noting that, in particular, Eq. ([T3|) describes 
Josephson-type vortices with size ^ A™ ^ I and flux 00- 

It follows from Eqs. pT|) and that equations de- 
scribing the coarse-grained phase ip{x) are the same as 
for Josephson junctions of conventional superconductors 
in the Meissner state. 

Specific symmetry of jc{x) might lead to a Fourier se- 
ries with some of the transforms being zero. In this case 
only the second term at the RHS of Eq. (fTT|) is non zero. 
In particular, if jc{x) is a stepwise function (see Fig. [1]), 
then we have i2k = and i2k+i = —2i/n (2k+l), where k 
is an integer. In the FS states with the field Ha located 
in the intervals AHa centered at "odd" resonant fields 
Ho]^ = (2k + l)Hi, Eq. ^ takes the form 



sinV', Aok = V7r(fc + l/2)Aj 



Similarly, for 
(fTTjl yields 



'even" resonant fields iJ( 



ck 



(15) 
2kHi, Eq. 



A^j^T/^" =sin2V', Aek 
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where (3k is a constant (/3i = 0.9, P2 = 1-1, = 0.9, and 
« 1 for fc > 3). It is seen from Eqs. ^ and HID 
that the "even" -field vortices carrying flux (j)o/2 are by 
the factor Xj/l ^ 1 wider than the "odd" -field vortices 
carrying flux (j)o- 

Next, we calculate the width of the intervals AHa of 
existence of the FS states by minimizing free energy, J-', of 
the junction. This can be done explicitly if we specify the 
spatial distribution of the critical current density, jc{x). 
Here, for brevity, we use the one harmonic model, i.e., 
we assume that jc{x) = ji sin(27ra;/?): 
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Hsn 
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where J- is normalized by H^^L/Sn, with Hsm = 
0o/47rAAm <C Hjn. If there is no vortices in the junc- 
tion, then we have ip — and 'p{x) — 2'K(j)iX / (pol + ipo, 
where (po is a constant. In this case the minimum value 
of the free energy, J-'m , is achieved for ipo — 6„i , 



27rA„ 
14>0 



sin [KN{m - (pi/cpo)] 
TTN{m - (pi/cpo) 



(18) 



where (j>a — 2XlHa. Finally, we fix the applied field Ha 
and minimize J-m with respect to the internal flux (pi. 
It follows from this calculation that if Ha is from the 
narrow interval uiHi — Hsm < Ha < mHi + ifsm, then 
the flux (pi is constant and equal to m0o, which means 
that Hi K, mHi. 

The above theoretical analysis can be supported by nu- 
merical simulations. To treat both statics and dynamics 




FIG. 3: (a.l)-(a.4) Magnetization, (j>i{(j>a), of a junction 
witfi stepwise critical current density and A'^ — 40; (a.l) 
the plateaus with (pi — const reveal the first four FS states. 
The curves (j>i[(j)a) in the vicinity of three FS states: (a.2) 
Ha = (a.3) Ha = 3Hi; (a.4) Ha = 5Hi. (b.l)-(b.2) 
Magnetization of a junction of conventional superconductors. 
(b.2) (j)i{(j)a) of a conventional junction in the Meissner state. 
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4: (a) The phase %l){x) distribution for vortices carrying 
■f>o at three "odd" resonant fields H\, 2>H\, and Siifi. (b) 



The same three curves are shown to collapse into one curve if 
the coordinate is normalized by Ai, A3, and A5. 
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FIG. 5: The phase '(/'(x) distribution for the first "odd", Ha ~ 
Hi, and "even" Ha = 2Hi resonant fields describing vortices 
with fluxes (^o {Ha = Hi) and (f>o/2 {Ha = 2Hi). 



of the phase difference f{x,t) we introduce time depen- 
dence into Eq. ([5]) and arrive to 



(p + aip — y?j^p" + icix) sin tp ~ {) , 



(19) 



where a ~ 1 is a decay constant. The second term in 
the LHS of Eq. (fT9|) describes dissipation. As a result of 
this dissipation the system ends up in one of the stable 
stationary states which is a solution of Eq. ^ . 

We used the finite-difference explicit method (see for 
details [l^) to solve Eq. (fT9|) assuming that the alternat- 
ing critical current density jc{x) is stepwise. The mag- 
netization curves (dependencies of 4>i on (j)a) obtained by 
numerical simulations are shown in Fig. [3l The plateaus 
with (pi = const are clearly seen for the series of the four 
FS states. The field inside the junction is constant for 
each of the FS states in contrast to junctions of con- 
ventional superconductors for which the internal field is 
constant only in the Meissner state. It is seen in Fig. 
El that the magnetization curve of a junction in the FS 
state is the same as the magnetization curve of a junc- 
tion of conventional superconductors but with the field 
Ha biased by vnHi. Our numerical simulations confirm 
that the width of the plateau in which the internal field is 
constant is proportional to i/jiml = Aj/Am as it follows 
from the above theoretical analysis. 

Next, we simulated numerically the dynamics and stat- 
ics of single Josephson-type vortices in the FS states (see 
Figs. mandO. We observed stable Josephson-type vor- 
tices carrying fluxes (/iq (for "odd" resonant fields) and 
(^o/2 (for "even" resonant fields). The width of vortices 
carrying flux 0o scales proportionally to -^/liml as follows 
from Eq. p5|) . Eq. ((T6)) shows that the width of vortices 
carrying flux 0o/2 is by factor Xj/l ^ 1 larger than that 

of 00- 

We found also that periodical motion of single vortices 
between the edges of the junction [x = and x — L) 
produce steps in the IV characteristics at a series of volt- 
ages Vm = 'm(j)QCs/L, where m is the number of fluxons 
inside the junction and Cg is the Swihart velocity [l6| . 
These high field steps are similar to the zero field steps 
in junctions of conventional superconductors [Til . 

To summarize, we find new equilibrium FS states in 



Josephson junctions with periodically alternating critical 
current density. The FS states exist if the applied field 
is from narrow intervals centered at equidistant series 
of resonant fields. If a junction is in FS state, the flux 
in the junctions is fixed and the maximum supercurrent 
across the junction is significantly enhanced. Two types 
of single high- field vortices with flux ipo or / and (/)o /2 exist 
in FS states. Periodic motion of these vortices between 
the edges of the junction results in high field steps in the 
IV characteristics. 

One of the authors (RGM) is grateful to J. R. Clem, V. 
G. Kogan, J. Mannhart, and C. W. Schneider for support 
and stimulating discussions. 
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